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Abstract 



^ . We calculate direct CP-violating observables in charged B — > VV decays arising from 



the interference of amplitudes with different strong and CKM phases. The perturbative 
strong phases develop at order as from absorptive parts of one-loop matrix elements 
of the next-to-leading logarithm corrected effective Hamiltonian. CPT constraints are 
maintained. Based on this model, we find that partial rate asymmetries between charge 
conjugate decays can be as high as 15-30% for certain channels with branching 
ratios in the 10~^ range. The small values of the coefficients of angular correlations, 
which we calculated previously to be of order 10~^, are not significantly degraded by 
the strong phases. The charge asymmetries of rates and angular distributions would 
provide unambiguous evidence for direct CP violation. 
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1 Introduction 



So far CP violation has been detected only in processes related to — mixing p| 
but considerable efforts are being made to investigate it in B decays. While the most 
promising proposal for observing CP violation in the B system involves the mixing 
between neutral B mesons 0, the particular interest in decays of charged B mesons lies 
in their possibilities for establishing the detailed nature of CP violation. Since charged 
B mesons can not mix, a measurement of a CP violating observable in these decays 
would be a clear sign for direct CP violation, which has been searched for in K decays 
with indefinite success as the measurements of e'/e do not yet exclude a zero value 0]. 

In non-leptonic charged B decays two main categories of CP-violating observables 
can be investigated: First, rate asymmetries P, |^, 

r - r 

«cp = Yqrf , (1) 

where F and f are the (partial) rates of the decay and its charge (C) conjugate, and 
second, azimuthal angular correlations J^, ||. 

The rate asymmetries occur even for spinless final states and require both weak and 
strong phase differences in interfering amplitudes. The weak phase differences arise 
from the superposition of amplitudes from various penguin diagrams and the usual W- 
exchange (if contributing). The strong phase is generated by final state interactions. 
At the quark level these strong interaction effects can be modeled by absorptive parts 
of perturbative penguin diagrams (hard final state interactions) |^ while predictions 
at the non-perturbative hadronic level are of course extremely difficult (soft final state 
interactions). Clearly we can not exclude that the weak transition matrix elements 
receive phases originating from soft final state interactions (resonances) between the 
produced vector particles. However, since the mass of the B is far above the usual 
resonance region, we expect these phase shifts to be small. 

It is clear that a significant contribution of penguin diagrams, and hence of the CKM 
IP] phase differences, is an exceptional case and requires either the absence or a strong 
CKM suppression of the tree contributions (as e.g. in charmless b s transitions). 

The rate asymmetries for exclusive two-body decays into pseudoscalars have been 
estimated by several authors using either the model of Bauer, Stech and Wirbel flOl 
(BSW) based on wave functions in the infinite momentum frame, or the perturbative 



methods developed by Brodsky et. al. The rate asymmetries acp can be quite 



large (of the order acp ~ 0.1) for some of the final states. However, the corresponding 
branching fractions of these decays are quite small, ranging from 10~^ (estimates with 
the BSW model [|l^) to 10~^ (estimates with the Brodsky-Lepage model |Tl|, [l^, |^). 
The magnitude of BR is therefore of the order of 10 ^ — 10 ^. 



2 



The second category of CP-violating observables involves the decay of the B meson 
into two vector particles B V1V2 with subsequent decays of Vi and V2 |0, §]• (In 
the following B will always denote the B~ meson and B its antiparticle.) By analyzing 
the azimuthal dependence of the vector meson decay products one can then isolate CP 
odd quantities. The advantage of this method is that the CP violating terms occur even 
when there are no strong phase differences between the interfering weak amplitudes. On 
the other hand these coefficients in the azimuthal correlation are also present when the 
CP-violating weak phase differences vanish. By measuring these coefficients in charge 
conjugate B^ decays one has the possibility to disentangle the effects of strong and 
weak phases [||. 

So far the "CP-odd" azimuthal angular coefficients have been calculated under the 
assumption that the strong phase differences vanish . These results give us an estimate 
of the effect we could expect from the CP-violating phase. Since strong phase differences 
are present, at least due to the hard final state interactions which lead to the rate 
asymmetries through the Bander-Silverman-Soni mechanism |p, we should include these 
strong phases also in the calculation of the "CP-odd" angular coefficients. 

The two categories for detecting direct CP violation are therefore complementary. 
The rate asymmetries occur only when non- vanishing strong phase differences are present 
but need no angular information. The azimuthal angular correlation terms need no 
strong phase differences but require joint angular distribution measurements of the de- 
cay products of the B mesons [14|. 

Although the explicit form of these angular distributions depends on the spins of 
the decay products of the decaying vector mesons Vi and V2, two formulas are sufficient 
to describe a general B— >VV decay p|. For instance, the angular distribution for the 
cascade decay B^ K*^p^ [K-k) (tt+tt^) has the following form: 

'^^^ .2/1 • 2/) , 2n 2q 

■ sm Ui stn 62 + — ■ cos Ui cos 62 



dcos9idcos92d(j) 4 F F 

-|-- sin26'i sin292[cii ■ coscp — (3i ■ sincj)] (2) 

+ -sin^ 9isin'^ 92 [c^2 ■ cos2(p — P2 ■ sin2(j) ] . 

In eq. 9i is the polar angle of the K in the rest system of the K* with respect to 
the helicity axis. Similarly 6*2 and (f) are the polar and azimuthal angle of the 7r+ in the 
p° rest system with respect to the helicity axis of the i.e. is the angle between the 
planes of the two decays K*~ Kit and p° tt+tt^. 
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The decay distribution is parameterized by the coefficients: 













r 




|/f+lP + |i/_iP 
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\Hq\^ + 


|//+lP + |/f_lP 


~ |Ho|2 + |//+l|2 + |J/_iP 



(3) 



where Hx = {Vi{\)V2{X)\TCeB\B) are the hehcity amphtudes (A = 0, ±1). Clearly eq. @ 
can be used also for all other decays where Vi and V2 decay into two pseudo scalar 
mesons. The other decay distribution, e.g. for B D*D* with D* Dg'y and 
D* — > Dtt can be found in ref. 0, where further examples are discussed. 

In our previous work, we have calculated the six angular coefficients and the branch- 
ing ratios for 36 decays B V1V2 with neutral and charged B mesons. Penguin con- 
tributions were taken into account only through leading-log short distance QCD effects 
and all strong phases were neglected. Non- vanishing "CP-odd" azimuthal angular corre- 
lations, i.e. coefficients /3i and P2 occurred only in the decays — * K*~uj, K*~p^, ujp~ 
(and the corresponding decays of In this work, we shall present results including 
strong phases from penguin diagram contributions to the matrix elements. We consider 
the rate asymmetries acp and the decay parameters from eq. for B^ and B^ decays, 
and investigate the influence of the strong phases on [3i^2- In addition, we base our treat- 
ment on the next-to-leading logarithmic short distance corrections evaluated by Buras 
et al. |]T3|, which is mandatory if one wants to systematically take into account the 
complete 0{as) penguin matrix elements. We include also some pure penguin modes 
which are of interest for the detection of CP effects via the rate asymmetry acp, and 
we give estimates of their branching ratios. 

The remainder of this paper is organized as follows. In section 2 we describe the 
effective weak Hamiltonian and the evaluation of the hadronic matrix elements. The 
CP- violating observables are discussed in section 3. The final results for the angular 
correlations and rate differences are discussed in section 4. Formulae for the matrix 
elements and some technical details about CPT cancellations can be found in the ap- 
pendices. 



2 The effective Hamiltonian 

2.1 Short distance QCD corrections 

For calculations of CP-violating observables it is most convenient to split the effective 
weak Hamiltonian into two pieces, one proportional to Vu = VubV*^ (or VubV*^ in the 
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case of b ^ d transitions) and the other one proportional to Vc = VcbV*^ (or VcbV*^ 
correspondingly) , 

(4) 



The two terms (g = u, c) 



i 

differ only by the quark content of the local operators, and for our purposes it is sufficient 



to consider only the following four-quark operators [ITB] : 



0[ 



(</) 



q' 



(5) 



q' q' 

where L and R are the left- and right-handed projection operators. The operators 
O3, . . . ,0q arise from (QCD) penguin diagrams which enter at order in the initial 
values of the coefficients, 



w) 



l+0(«,) (« = 2) 
0{as) (otherwise) 



or through operator mixing during the renormalization group summation of short dis- 
tance QCD corrections. The renormalization group evolution from fi ^ Mw to yU ~ rrib 
has been evaluated in next-to-leading logarithmic (NLL) precision by Buras et al. |jl5 



These authors also demonstrated how the 0{as) renormalization scheme dependence 
can be isolated in terms of a matrix r^j by writing 



Si 



Att ''''' 



(6) 



where the coefficients c,- are scheme independent at this order. The numerical values 



for A 



MS 



350 MeV g, = 150 GeV and fi = rrib = 4.8GeV are 



-0.324 
0.017 , 
0.011 , 

Contributions from the color magnetic moment operator 



Cl 

C3 
C5 



C2 

C4 

C6 



1.151 , 
-0.038 , 
-0.047 . 



(7) 



On 



9s 



167r2 



^ This value of A^-L- translates to w 250MeV, which is about the value from a recent compi- 

MS MS ' ^ 

lation of G. AltareUi (a1£L = 240 ± 90MeV |l|]) 
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with a coefficient of the order —0.15, will allways be neglected in the following, because 
already its tree level matrix elements are suppressed by a factor as/47r and it cannot 
provide interesting absorptive parts in the decays considered here. 



2.2 Quark-level matrix elements 

Working at NLL precision, it is consistent - and necessary in order to cancel the scheme 
dependence from the renormalization group evolution - to treat the matrix elements of 
Ties at the one-loop level. These one-loop matrix elements can be rewritten in terms of 
the tree-level matrix elements of the effective operators, and one obtains: 



{sq'q'\niS\b) = Y.c^{^^) 



The functions niij are determined by the corresponding renormalized one-loop diagrams 
and depend in general on the scale /i, on the quark masses and momenta, and on the 
renormalization scheme. The various one-loop diagrams can be grouped into two classes: 
vertex- corrections, where a gluon connects two of the outgoing quark lines (fig. la), and 
penguin diagrams, where a quark- antiquark line closes a loop and emits a gluon, which 
itself decays finally into a quark-antiquark pair (fig. lb). 

When expressing the rhs of eq. (|^) in terms of the renormalization scheme indepen- 
dent coefficients q, the effective coefficients multiplying the matrix elements {sq'q'\Oj^^\by^^'^ 
become 

a 



c'^ ^c- 



+ (mii-rij) . (9) 



The renormalization scheme dependence, which is present in niij and rjj, explicitly 
cancels in the combination rrijj — Tij. This reflects the familiar fact that one- loop matrix 
elements have to be included in order to compensate for the order a.^ renormalization 
scheme dependence which enters through the coefficients q in Hes (generated by the 
NLL renormalization group evolution). For instance, the scheme dependence in the 
coefficients of the penguin operators, which enter via Cj{sq'q'\Oj'^^\by'^^^ for j = 3, 4, 5, 6 
in eq. (H), is cancelled by penguin-like matrix elements (fig. lb) of the operators. 

The effective coefficients cl% receive contributions only from vertex-correction dia- 
grams, which will not be included in the following (see the discussion in Section 3). For 
a general SU{N) color group the remaining effective coefficients can be brought into 
the following form 



1 a.s 



of = ^"3-^^(Q + c,) + 

Cf = C4 + - — (Q + Cp) + ••• 
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Cf = C5 



1 a. 



{ct + Cp) + 



„eff 



2iV47r 



(10) 



where we have separated the contributions Ct and Cp from the "tree" operators Oi^2 and 
from the penguin operators O3...6, respectively. The eUipses denote further contributions 
from vertex- correction diagrams. 

In addition to the contributions from penguin diagrams with insertions of the tree 
operators Oi'l 



Ct = C2- 



10 2 ml , ^ / A;2 X 



9 J 



where AFi is defined in appendix A, we have evaluated the penguin diagrams for the 
matrix elements of the penguin operators: 



C3 



280 2 ml 2 ml fk\ fk'^\ 



+ (c4 + ce) • Yl 

j=u,d,s 



11^ 3 
"9 



'■mi 



<mt 



Q 3 //^ ^m-^ 



(12) 



2.3 Hadronic matrix elements in the BSW model 



In order to evaluate the hadronic matrix elements of TYefr we represent the helicity 
amplitudes in terms of three invariant amplitudes, a, h and c: 

h 



= ei^(A)*e2.(A)^ 



ic 



(13) 



mim2 mim2 

where pi^2 and mi^2 are the four- momenta and masses of Vi,2, respectively. 

The coefficients a, h and c have strong phases 5 from final state interactions (e.g. 
between the two vector particles Vi and V2) and weak phases originating from the 
CP violating phase in the CKM matrix. In general, the invariant amplitudes are a 
sum of several interfering amplitudes, a^, 6^, and c^, respectively, corresponding to, for 
instance, various different isospin contributions. Then the phase structure of a, h and c 
is: 



(14) 



The helicity amplitudes Hx for the decay of i? — > V1V2, where and V2 are the 
antiparticles of Vi and V2 respectively have the same decomposition as (12) with a — > 
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a, 6 — i> 6 and c —c. The amplitudes a, b and c have the analogous phase structure as 
above with the weak phases changing sign (pak,bk,Ck —(pak,bk,CkJ i-^- ^k, bk and Ck 
have the same strong phase shifts and the opposite weak phase compared to a^, bk and 

Ck- 

From the decomposition eq. ([131) one finds the following relations between the helicity 
amplitudes and the invariant amplitudes, a, b, c: 

H±i = a ± cV — 1 and Hq = —ax — b (^x'^ — 1^ , (15) 

where 

P1P2 '^R ~ ''^1 ~ ""^i 



771x1712 2mim2 

When there are no strong interaction phases, one has a = a*, b = b* and c = c*. Due 
to the sign change in front of c in Hx, we have in this case 

fj TJ* fj TT* 

U±l — , — . 

To take into account long distance QCD effects which build up the hadronic final 
states, we follow Bauer, Stech and Wirbel [1^: With the help of the factorization hy- 



pothesis ||18[ the three-hadron matrix elements are split into vacuum-meson and meson- 
meson matrix elements of the quark currents entering in Oi, . . . , Og. In addition, OZI 
suppressed form factors and annihilation terms are neglected. In the BSW model, the 
meson-meson matrix elements of the currents are evaluated by overlap integrals of the 
corresponding wave functions and the dependence on the momentum transfer (which 
is equal to the mass of the factorized meson) is modeled by a single-pole ansatz. As 
a first approximation, this calculational scheme provides a reasonable method for esti- 
mating the relative size and phase of the tree and penguin terms that give rise to the 
CP-violating signals. 

Concerning the QCD coefficients and how 1/A^ terms are treated, it is well known 
TP| that this model has problems accounting for the decays with branching ratios which 



are proportional to where 

1 _ 

02 = Ci + — C2 . 

This is due to the fact that 02 has a rather small value |a2| ~ 0.06 when using the short- 
distance QCD corrected coefficients. An analogous effect is also known in nonleptonic D 
decays fl^, and several authors advocated a modified procedure to evaluate the factor- 



ized amplitudes |]T0|, There, only terms which are dominant in the 1/A^ expansion are 
taken into account. Recently there has been much discussion in the literature concern- 
ing these issues. Some authors have argued that QCD sum rules validate this procedure 



21| . As our model for evaluating the matrix elements of the weak Hamiltonian we also 
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choose this leading approximation and use the QCD corrected coefficient functions 
Cj given above. We note that the terms proportional to in eq. ([T0|) must then be 
dropped as well. 

The strong phase shifts are generated in our model only by the absorptive parts (hard 
final state interactions) of the quark-level matrix elements of the effective Hamiltonian. 
Of course, when factorizing the hadronic matrix elements, all information on the crucial 
value of the momentum transfer k"^ of the gluon in the penguin diagram (fig. lb) is lost. 
While it has been attempted |]12| to model a more realistic momentum distribution by 



taking into account the exchange of a hard gluon, we will use here for simplicity only a 
fixed value of k"^. From simple two body kinematics [El] or from the investigations in 



ref. [12] one expects k to be typically in the range 



^ ^ P ^ ^ . (16) 
4 2 ^ ^ 



3 CP-violating observables 

In general, CP-violating observables require the interference of two amplitudes with 
different weak phases, 0, from the CKM factors. To investigate the necessity of strong 
phases for the observables in — > VV decays (see section 1) and their interplay with 
the weak phase factors, we decompose the amplitudes into contributions proportional 

to Vu and Vc 

/l = ^„.AH+y^.AW , (17) 

where A stands for a generic decay or helicity amplitude. The amplitudes for the CP 
conjugate process, are then obtained from eq. (0) by replacing Vq by v* 

Observables which involve only the real parts of the interfering amplitudes, like the 
decay rate or the parameters ai^2, can signal CP violation only when one compares 
them with the corresponding quantities of the charge conjugate decay channel, and 
when both, non-vanishing weak phase differences and strong phase shifts, Su — Sc, are 
present. For instance, (defining Vg = |fg|e*'^'') 

F - r ~ Im[vuV*] ■ Im[AuAl] ~ sin(0„ - 0^) ■ sin((5„ - 4) . (18) 

On the other hand, the decay parameters (3i {i = 1,2) can have non-zero values 
in the presence of either weak or strong phases alone. Then, by comparison with the 
parameters Pi (/3f ^) of the C (CP) conjugate decay, one can, in principle, establish a 
weak phase difference even for vanishing strong phases 

A + A = f^i-f^?^ ~ Im[vuV*] ■ Re[AuAl] ~ sin(0„ - 0^ ■ cos(5„ - 4) . (19) 
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or measure the strong phase shifts even for neghgible weak phases 



A - A = A + A ~ M^A] ■ iHAjAl] ~ cos(0, - ■ sin(5,- - 6k) , (20) 

where Vj and Vk are not necessarily different. In eqs. (|l^) and (^) we have dropped 
terms proportional to Im[vuV*] ■ Im[AuA*^ which arise from the different denominators, 
r and f , in A and A- Note also the relative sign between A and due to the parity 
reflection. 

If there are no CP- violating weak phases then A = ~A; Q^i = "Si, and F = F, while 
the absence of strong phases implies A = A (and, of course, = cxi and F = f). 
Interesting CP differences in the case B VV, which do not require strong phases 
[see eq. (|T^)] and which are proportional to weak phase differences, are [using the phase 
definitions of eq. (p^j 

Im{H+iH*_^ + H+iH*_^) = -4Va;2 - l^l s^n{4>ai - <Pcj) cos{5ai - 6cj)\aiCj\ (21) 

and 

Im{H+iH* - H^iH* + H+^H* - H^^H*) = 

-4(x^ - l)t ^ sin{(j)ci - 4>bj) cos{6ci - Sbj)\cibj\ (22) 



-AxVx"^ -lYl Sin{(j)ci - (paj) C0s{5ci - 5aj)\Ciaj\ . 

hj 

Here, are the amplitudes of the charge conjugated process and the subscripts of the 
weak and strong phases refer to different weak (or isospin, etc.) contributions. 

The presence of strong phases is unambiguously demonstrated by a partial rate 
asymmetry as well as angular correlations of the following kind: 

- = -("2 - a2)cos2(j) - (A - A)sm20 (23) 

L a(p 1 (30 

Other terms can be isolated by examining the dependence of the differential rate 
difference between same hemisphere (SH) events (e.g. < 6*1,^2 < f) and opposite 
hemisphere (OH) events (e.g. < ^1 < |, | < ^2 < tt): 

In general the dominant terms in the angular correlations are F^/F, F^/F, ai and 02. 
The terms A and A are small and they are nonvanishing only if the helicity amplitudes 
H^i,H_i and Hq or the invariant amplitudes a,b and c, respectively, have different 
phases. In all those channels where factorization is possible only in one way, this is not 
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the case because all matrix elements become simply proportional to each other; this is, 
of course, due to the simplicity of our model. For the same reason one has = cij (and 
F-r/r = Ft/F) in all the channels with vanishing (3i: The overall (weak and/or strong) 
phase factor cancels in the ratios that enter in the definition of the [see (^]. 

While the next-to-leading logarithmic precision of the effective Hamiltonian allows 
one to consistently calculate all amplitudes at order as and to include all one-loop matrix 
elements, some care is necessary when evaluating CP-violating asymmetries of the decay 
rates or of the observables of the angular distribution. In particular, one should make 
sure that the rate asymmetries for sufficiently inclusive channels remain consistent with 
CPT constraints in certain mass limits [Q. 

In order to specify a procedure which meets this requirement, we recall that absorp- 
tive parts, which arise from intermediate states having the same quark-gluon content as 
the final state, play a crucial role for CPT consistency of rate asymmetries in inclusive 



decays: It has been shown (see e.g. ^) that the inclusive rate difference due to ab- 
sorptive parts from these "fiavour-diagonal" interactions must, and in fact does, cancel 
if all interferences are taken into account which contribute at a given order in (see 
appendix B for an example of such a "CPT cancellation"). Of course, this cancellation 
need not be generally true when the phase space of the final state is restricted or when 
exclusive decays are considered. Nevertheless, in our model for the exclusive amplitudes 
these CPT cancellations hold for many of the diagrams due to the factorization of the 
hadronic matrix elements. We shall assume that analogous CPT cancellations are — 
at least approximately — valid for all interferences of exclusive amplitudes to be eval- 
uated here, and we therefore neglect absorptive parts from fiavor diagonal rescattering 
throughout. 

In this approximation all imaginary parts of the terms AFi{k'^ /m?^ are to be dropped 
in (|lTD and ([l^) when q refers to the fiavor of a gg-pair which is present in the final 
state. Moreover, no absorptive parts of vertex correction diagrams have to be evaluated, 
because they are always fiavour diagonal. For consistency, the same procedure as for the 
rates should, of course, be systematically applied when evaluating the decay parameters 
(§) and their CP-differences. To get an idea of the quality of our approximation, we 
calculated the flavour-diagonal absorptive parts for the case of all penguin-like matrix 
elements, (0) and ([T2|) , and we found that their effect is indeed small. 

In our calculation we do not explicitly drop higher order terms which arise, for 
instance, through interferences among (real and imaginary parts of) the order as matrix 
elements. However, such terms can not introduce the above mentioned inconsistencies 
with CPT because the flavour-diagonal absorptive parts are discarded. A completely 
systematic treatment of the higher order terms, some details of which we describe in 
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appendix C, would require to expand all products of interfering amplitudes — and not 
the amplitudes themselves — in term of and the other couplings of the (effective) 
theory. 



4 Results and Discussion 

For a numerical analysis of the decay parameters and their CP-violating effects within 
our model, we need to specify the CKM matrix elements and the current form factors. 
It is well known that fits for the parameters f\ 

p = COS(5i3 Si3/(Si2S23) 

T] = sinSi3 Sis/{si2S23) (25) 

of the CKM matrix depend critically on the value of the 5-meson decay constant Jb- 
The solution for lower fs values leads to a negative p while higher fs values render p 
positive. 

We have calculated our results for the positive p solution, with the values 

p = 0.32, 7] = 0.31 



(i.e. Si3 = 0.0045, 613 = 44deg) from the analysis by Schmidtler and Schubert |p5|| for 
Jb = 250 MeV (giving rrit = 135 ± 27 GeV). A more recent analysis by Ali and London 
p7| based on the latest information on Vub yields similar results. For comparison, we 
will also show the asymmetries calculated with the negative p solution: p = —0.41 and 
T] = 0.18 corresponding to fs = 125 MeV (giving rrit = 172 ± 15 GeV). 

The main purpose of this work is to calculate the effect of (perturbative) strong 
phases on the angular correlation coefficients obtained earlier with the strong phase 
put to zero Moreover, we have included pure penguin channels and calculated 
partial rate asymmetries. (The pure penguin modes are p~K*^, K*^^, and 

—>■ K*^K*^.) Of course, the estimates of the CP-violating observables given here 
may suffer from large uncertainties due to strong phases from soft final state interactions 
and, therefore, can at most be indicative of the expected orders of magnitude. 

First we consider the results for the p positive case (tables 1 and 2). To see the 
effect of strong phases we also show in the following tables the results without strong 
phases from imaginary parts of the one-loop matrix elements (the values are given in 
parentheses and only where different from those with strong phases included). These 
numbers may differ from the results in M, R] since we include here throughout the order 



These coincide with the parameters p and rj of the Wolfenstein representation for small angles 
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as real parts of the penguin-like one-loop matrix elements. The branching ratios are 
calculated using tb = 1.49 psec, and a dash ( — ) in the tables indicates values which 
must be exactly zero in our model. 

The parameters /3i and (32 of the azimuthal decay distribution are only non-zero for 
the decays B~ —>■ K*~uj, K*~ and p~uj. The other decays have = in our model 
because all three helicity amplitudes have the same overall (combined weak and strong) 
phase. This is always the case when there is only one way to factorize the hadronic 
matrix element, e.g. in K*^J/'^, D^^D*^ , etc., and in all pure penguin modes. 

A special case is B~ P^P^'- if isospin breaking (due to the mass difference between p~ 
and p^) is neglected, no penguin contributions are present^ and, hence, no weak phase 
differences can occur. Concerning the CP-violating effects in the angular distributions 
the most promising decays in the case of vanishing strong phases (values in parentheses) 
are B~ K*~uj and B~ — > K*~ p^: They have branching ratios of the order (1-3) x 10~^ 
and \(3i\ ^(l-4)xlO~^ whereas the (32 are a factor of ten smaller. (See table 1.) 

Table 1 also illustrates the influence of the strong phases generated by the absorptive 
parts of the matrix elements. They have two effects: First, they produce the rate 
asymmetry acp 7^ 0, which is given in the third column, and second, they generate 
different angular distributions for certain charge-conjugate decays, i.e. 7^ and 
(3i + A- 

The rate asymmetry acp is appreciable in some of the cases, e.g. for B^ —>■ K*^uj 
(acp ~ 28%) and B~ K*" p^ (acp ~ 15%). Both decays have approximately the same 
branching ratio of the order (2-4) x 10"^. For the pure penguin modes these asymmetries 
are either smaller (^ 1%) or, as in the case of B^ K*~K*^, the branching ratios are 
tiny (~ 10^^). Interesting is the decay B^ D*^D*^ with a branching ratio of 0.1% 
and a rate asymmetry^of about 1%. For B^ — * p'K*^ the branching ratio and acp are 
similar to the results for tt~K*^ obtained recently by Fleischer 



For (3i and /?2 "we observe a significant effect from the strong phase shifts for the most 
interesting case B^ K*~uj. From our results for (3i and (32 for the charge conjugate 
decay B^ K*^uj (see the corresponding line in table 1) we find that (3i + /3j is not 
drastically changed as compared to the case with no strong phases. We obtained from 
table 1 (/3i + (3i)/2 = -29 x 10"^ and {(32 + (32) /'^ = 2.8 x lO^^ which is to be compared 
to = (3i = -37 X 10-^ and p2 = (32 = 3.6 x 10"^ On the other hand, we find 
iPi - (3i)/2 = -14 X 10-3 and {p2 - (32)/2 = 1.4 x 10"=^ which shows the effect of 
the strong phases [see eq. (16) and eq.(17)]. In the case of K*~p^ the behavior is less 
significant, e.g. {(3i + (3i)/2 = 6.7 x 10-^ and {(3i - (3i)/2 = -0.4 x 10^^ (compared to 

^This has been missed in 

^ This channel is rather insensitive to the particular choice of fc^ because the leading absorptive part 
comes from a uu-cut with a low threshold. 
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A = /?! = 7.2 X 10 without strong phases). In the puj channel the pattern is even 
more pronounced: while the (3i + (3i are again not changed substantially by the strong 
phases, the differences, -/3i)/2 = 0.16 x lO^^ and (/32 -/32)/2 = -0.002 x 10-^ are 
now larger than the sums. 

In the analogous calculations for ai — ai (and r^r/r — Tt/'^) we find significant C- 
differences only in the case of K*~uj, where the relative differences of Tt/T, ai and a2 
are of the order of several percent (detailed numbers can be extracted from table 1). In 
the two other cases, K*~p^ and p~uj, these differences are smaller. 

We have also performed the corresponding calculation with the 1/A^ terms included 
throughout (see table 2). As a result, some of the predicted branching ratios (BR) change 
drastically; for instance, BR{B~ —>■ K*~ J/'^) is decreased from 3.8 x 10~^ to 1.6 x lO"'^. 
Since the average experimental value of this branching ratio is (0.17 ± 0.05)% ||l9| the 
version without 1/A^ terms is to be preferred in this decay. The pure penguin modes 
of the h s transitions are rather insensitive to the inclusion of the 1/A^ terms with 
the exception of K*$ whose rate is increased by a factor of three. The branching ratios 
for all other b —>■ s transitions in the upper part of the tables remain more or less 
unchanged, while (3i for B~ K*~uj and B~ — > K*~p^ become smaller and look less 
interesting. The pattern for h ^ d transitions (see the lower part of table 2) is different 
when compared with table 1: the branching ratio for the pure penguin mode — 
is drastically reduced and also the rate for B^ —>■ p^J/'^ becomes more than an order 
of magnitude smaller. 

It is obvious that the color unsuppressed decays B^ D*^D*^, D*^D*^ (with 
amplitudes proportional to C2 + Ci/N) are not significantly influenced by the treatment 
of the 1/A^ terms. We should also note that the coefficients Tt/T, Tl/T, ai and a2, in 
the angular distribution are not very sensitive to the treatment of the 1/A^ terms and 
to assumptions about the strong and weak phases. They depend mainly on the helicity 
structure of the matrix elements and they are therefore more important for testing the 
underlying model assumptions, in particular, the current matrix elements in the BSW 
model in conjunction with the factorization hypothesis. 

To give an impression how sensitive our results are with respect to the solution am- 
biguity in the CKM parameter determination we have repeated the calculations with p 
negative. The results (without 1 /N terms) are shown in table 3 and should be compared 
with the numbers in table 1 (with strong phases). Generally the branching ratios and 
asymmetries are similar in magnitude. In the interesting K*~uj and K*~p^ final states 
the branching ratios decrease or increase, respectively, while the rate asymmetries and 
azimuthal asymmetries vary in the opposite direction. The pure penguin modes change 
less but in a similar way. Since the charge asymmetries of the various observables can 
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be more pronounced for a p negative CKM matrix, the interesting charge conjugate 
channels are also included in table 3. 

We mention that some of the pure penguin modes have been calculated by other au- 



thors. Davies et al. obtained a comparable result for the partial rate oi B ^ K*(^ |29 



using only leading logarithmic order QCD coefficients. Dong-sheng Du et al. |20| calcu- 
lated the p$ rate using also coefficients of Buras et al.. They obtained similar results, 
in particular also the strong suppression of this rate when 1/N terms are included. 

In this work we have confined our attention to direct CP violation in decays of 
charged B mesons. The strong phases also contribute to some neutral B decays, where 
their effect on the CP-violating time dependence is a further complication in the analysis 
of CP violation arising from the interference of mixing and decay amplitudes. This 
subject is currently under investigation. 

Finally, we would like to note that the partial rate asymmetries we have found are 
larger than those reported in B decays to two pseudoscalar mesons, making B — > VV 
an attractive channel for probes of direct CP violation. 
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Appendix A: Matrix elements 

The momentum dependence of the penguin-hke matrix elements entering in and 
(1^) is given by 

AFi(z) = -4 [\{l-u)in[l - zu{l-u)]du 
Jo 

where, setting r = — 4/ z\, 

'r-Cngl (;<0) 

—TTT + 2rarctan r {0 < z < 4) 

+irTT + rinY^ {z > 4) 



For completness, we hst here also the factorized matrix elements of Ties for the pure 
penguin modes not yet presented in former publications 0: 

{p-K*'\n^^\B-) = v,V2Gp{j^cf + cf){p'\d^,bL\B-){K*'\-sYd\0) , (27) 

{K*-m^S\B-) = v,V2G^{{l + j^)cf + {l + ^)cf + cf + j^cf) 

■{K*-\s^^bL\B-){<!>\sYs\0) , (28) 
(p-m'eSlB-) = v,V2G^{cf + ^cf + cf + ^cf) 

■{p-\dj,bL\B-){<!>\srs\0) , (29) 
{K*-K*'\Hi^\B-) = v,V2Gp{j^cf + cf){K*-\s^,bL\B-){K*''\dYs\0) . (30) 

Appendix B: CPT cancellations in inclusive rate asymmetries 

In this appendix we shall illustrate a typical example of a "CPT cancellation" for the 
CP-violating rate difference in inclusive decays. In contrast to the discussion within 
the full theory of refs. ||^, |22[ we adopt here the framework of the effective Hamiltonian 
where one can perform a completely analogous analysis. To render the correspondence 
to diagrams in the full theory more obvious, we assume the simplified situation where 
all Ci are zero, except C2. 

An illustrating channel is the charmless Ab = As = 1 transition b — > sum At 
order as the rate difference acp ~ F — f can arise only from the interference between 
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the tree-level matrix element of and the penguin diagram of (see fig. 2). The 
resulting asymmetry is small because the absorptive part of the penguin is kinematically 
suppressed by the cc-threshold. At order a1 one finds two types of contributions: First, 
the interference between penguin- like matrix elements of Oq^^ and (see fig. 3a), and 
second, interferences between the tree level matrix element of 0^2^ and order g'^ diagrams 
with an insertion of 02^. In particular, there is the interference with a penguin-like 
diagram of 02^ having an additional vacuum polarization on the gluon line (fig. 3b). 
An absorptive part of this diagram can be generated through a ww-pair inside the loop. 
In this case, the two order a1 interferences depicted in fig. 3 differ only by interchanging 
the role of the final and the "cut" state, which generates the absorptive part. Applying 



the Cutkosky rules one readily finds (see also ||23| , |22| ) that both interferences have the 



same size, and their combined effect cancels due to the relative minus sign in 

r - f ~ Im[AuAl] = Im[Au] ■ Re[A,] - Re[Au] ■ Im[A,] . 

Note that the cancellation between the diagrams of, for instance, fig. 3 works also for 
the (factorized!) exclusive amplitudes because the quark loop contributes in both cases 
the same multiplicative factor. 



Appendix C: Systematic expansion of the observables 

In the framework of the effective Hamiltonian described in section 2, the natural choice 
for the couplings, in terms of which the observables are to be expanded, is at = 
together with the manifestly renormalization scheme independent coefficients q. Of 
course, CPT cancellations analogous to the ones discussed above are present for each 
interference proportional to CjCjO" with n > 2. In table 4 we have specified the order 
up to which our calculation of the various observables is complete. Here, we use q and 
Cp as generic notations for any of the coefficients {ci,C2} of Oi^2 or of the coefficients 
{c3, . . . , Cq} of the penguin operators, respectively [not necessarily in the combinations 
defined in (0) and ([ID]. 

For decay modes with tree and penguin contributions we do not cover any terms 
of order in a complete way. However, since we are counting the powers of ^^(mb) 
within the effective theory, the interference between local penguin contributions (~ Cp) 
and absorptive penguin-like matrix elements (~ a^Q) is included here already at order 
(while it is order in the full theory). For instance, the local parts of the subdiagrams in 
fig. 3 (viewed as diagrams in the full theory), which are enclosed by a box, are taken into 
account in the effective theory through the penguin operators O3, . . . ,0q. Although the 
coefficients Cp, being implicitly of order as{Mw), are small, the order a^QCp interferences 
can be numerically important for the rate asymmetries in some charmless decays where 



17 



the ascf contributions (see e.g. Fig. 2) are kinematically suppressed due to the cc- 
threshold 0. 

In the case of pure penguin modes, i.e. final states which have at tree level only 
contributions from the penguin operators O3, . . . , Og, interferences of order alcf can 
arise only as a product of two order g"^ amplitudes. Therefore, the penguin-like one- 
loop matrix elements are sufficient for the complete treatment of these interferences. 
On the other hand, we may neglect systematically all terms of order a^Cp and a^CfCp. 
Thus, we retain effectively the terms which originate from order contributions in the 
full theory (and before the renormalization group evolution which sums up powers of 
as X CniJ? / M^r) . Of course, this procedure is also numerically sensible, because q ^ Cp 
by at least a factor of five (~ l/a^) due to their origin from penguin diagrams. 

At any order higher than the ones listed in table 4, real parts of vertex-correction 
diagrams are necessary for a complete treatment. Their imaginary parts, which enter 
already in the terms enclosed by parentheses in the table, correspond to fiavour-diagonal 
rescattering and are neglected throughout our calculations. Of course, a complete treat- 
ment of all order interferences also requires one to calculate a large number of two-loop 



diagrams |^ for the matrix elements of Ties- In fact, the next-to-leading logarithmic 
precision for the coefficient functions, just allows one to evaluate their absorptive parts 
(but not their real parts) without encountering new renormalization scheme dependen- 
cies. 
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Figure captions 



Fig. 1: Two types of one-loop matrix elements: (a) Vertex corrections, and (b) penguin 
diagrams. The square box denotes an insertion of one of the four-quark operators 
ofeq. (D. 

Fig. 2: Two interfering diagrams generating a rate asymmetry at order a^. The ab- 
sorptive phase arises from the cut state indicated by the dashed line. 

Fig. 3: Examples of order interferences: (a) Two 0{gl) penguin diagrams, and (b) 
a tree and a 0{gf) diagram. The rate differences due to the absorptive parts from 
the uu cut (dashed line) cancel. The local part of the subdiagrams enclosed by 
the dotted box is taken into account by penguin operators in the effective theory. 

Table captions 

Tab. 1: Branching ratios, rate asymmetries and angular correlation coefficients, using 
matrix elements without 1/N Terms for a p positive CKM Matrix (p = 0.32,// = 
0.31). Values in parentheses correspond to the case without strong phases. 

Tab. 2: Branching ratios, rate asymmetries and angular correlation coefficients, using 
matrix elements with 1/N Terms for a p positive CKM Matrix (p = 0.32, t] = 0.31). 
Values in parentheses correspond to the case without strong phases. 

Tab. 3: Branching ratios, rate asymmetries and angular correlation coefficients, using 
matrix elements without 1/N Terms for a p negative CKM Matrix (p = —0.41, rj = 
0.18). 

Tab. 4: Orders of as, Ct and Cp at which the treatment of the various observables is 
complete (neglecting absorptive parts from flavour-diagonal rescattering in vertex- 
correction diagrams). 
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Table 1 



Matrix Elements without 1/N Terms and with (without) Strong Phases 




p positive CKM Matrix: p 


= 0.32,r? = 


0.31 






Channel 


BR 


acp 


Ft 

r 


oil 






(32 






[%] 








[10-3] 


[10-3] 




6 — > s transitions: 


Ar — n 

LAC — U 


/\h — Ac 


— 1 

— 1 








1.3x10"^ 


0.54 


0.107 


— U.oo4 


U.UUy 








(1.2x10-5) 


(-) 












r> — > K +9 


5.5x10-6 


1.22 


0.137 


— 0.080 


0.017 








(4.3x10-6) 


(-) 












R~ K*~ -V- ij] 


2.4x10-6 


28 


0.089 


-0.316 


01 1 


—44 


4 9 


B+ K*+ + oj 


1.3x10^6 




0.073 


-0.300 


0.012 


-15 


1.5 




(1.7x10-6) 


(-) 


(0.080) 


(-0.307) 


(0.012) 


(-37) 

V / 


(3.6) 


B' K*- + 


3.7x10-6 


15 


0.105 


-0.332 


0.009 


7.1 


-0.63 


B+ K*+ + pO 


2.7x10-6 




0.103 


-0.330 


0.009 


6.3 


-0.55 




(3.0x10-6) 


(-) 


(0.104) 


(-0.331) 


(0.009) 


(7.2) 


(-0.64) 


B- K*- + J/* 


3.8x10-3 




0.427 


-0.621 


0.123 








(3.8x10-3) 


(-) 












B^ D*~ + D*^ 


2.4x10^2 


-0.05 


0.476 


-0.664 


0.183 








(2.4x10-2) 


(-) 














6 — > d transitions: 


As = Ac = 0, Ab 


= 1 






£> — ' p T" Vf 


1.1x10-^ 




0.130 


-0.364 


U.Uli 








(l.lxlO"'^) 


(-) 












B- K*- + K*^ 


4.2x10-^^ 


-18 


0.112 


-0.352 


0.014 








(4.8x10-^) 


(-) 












B~ ^ p~ +L0 


1.3x10-5 


-6.6 


0.083 


-0.298 


0.007 


+0.22 


-0.003 


B+ ^ p+ +u 


1.5x10-5 




0.083 


-0.298 


0.007 


-0.10 


+0.001 




(1.4x10-5) 


(-) 








(0.05) 


(-0.001) 


B-^p-+ pO 


1.3x10-5 




0.083 


-0.298 


0.007 








(1.3x10-5) 


(-) 












B- ^ p- + J/^ 


1.8x10-^ 




0.388 


-0.597 


0.097 








(1.8x10-^) 


(-) 












B- D*° + D*- 


1.2x10-3 


0.89 


0.456 


-0.660 


0.172 








(1.2x10-3) 


(-) 
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Table 2 



Matrix Elements with 1/N Terms and with (without) Strong 


r lidoco 






p positive CKM Matrix: p - 


= 0.32,7? = 


0.31 






Channel 


BR 


acp 


r 


ai 












[%] 








[10 3] 


r-i r\ 

[10 3] 




6 — > s transitions: 


Ac = 0, A6 = As 


= 1 






B- ^ p- + K*^ 


1.0x10"^ 


0.56 


0.107 


-0.334 


0.009 


— 


— 




(9.3x10-6) 


(-) 












B~ K*- + $ 


1.5x10-^ 


0.56 


0.137 


-0.385 


0.017 


— 


— 




(1.4x10-^) 


(-) 












B~ K*- + u 


2.6x10-6 


29 


0.107 


-0.334 


0.009 


-0.95 


0.09 




(1.8x10-6) 


(-) 


(0.107) 


(-0.333) 


(0.009) 


(-1.2) 


(0.12) 


B- K*- + p° 


2.5x10^6 


30 


0.106 


-0.333 


0.009 


-1.6 


0.14 




(1.8x10-6) 


(-) 


(0.106) 


(-0.333) 


(0.009) 


(-1.9) 


(0.17) 


B- K*- + J/^- 


1.6x10""^ 




0.427 


-0.621 


0.123 


— 


— 




(1.6x10"'^) 


(-) 












B- D*~ + 


2.0x10-2 


-0.05 


0.476 


-0.664 


0.183 


— 


— 




(2.0x10-2) 


(-) 














h ^ d transitions: 


As = Ac = 0, A6 


= 1 






B- ^ p- + ^ 


1.6x10-11 




0.130 


-0.364 


0.011 


— 


— 




(1.6x10-11) 


(-) 












B- K*- + K*^ 


3.1x10-^ 


-18 


0.112 


-0.352 


0.014 








(3.6x10-^) 


(-) 












B- ^ p- +UJ 


2.5x10-^ 


-4.1 


0.084 


-0.299 


0.007 


0.23 


-0.003 




(2.6x10-^) 


(-) 








(0.17) 


(-0.002) 


B- ^ p- + 


2.3x10"^ 




0.083 


-0.298 


0.007 








(2.3x10-5) 


(-) 












B' ^ p' + J/* 


7.0x10-6 




0.388 


-0.597 


0.097 








(7.0x10-6) 


(-) 












B- D*^ + D*- 


1.0x10-3 


0.88 


0.456 


-0.660 


0.172 








(1.0x10-3) 


(-) 
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Table 3 



Matrix Elements without 1/N Terms and with Strong Phases 




p negative CKM Matrix: p = 


= -0.41,7] 


= 0.18 






Channel 


BR 


acp 


r 


a\ 






(h 






[%] 








[J-U J 


[10-3] 




b ^ s transitions: 


Ac = 0, A6 = As = 1 






B- ^ p- + K*^ 


1.3x10"^ 


0.33 


0.107 


-0.334 


0.009 


— 




B- K*- + $ 


5.2x10-6 


0.75 


0.137 


-0.385 


0.017 


— 




B- K*- +u 


6.5x10-^ 


83 


0.135 


-0.342 


+0.004 


-131 


—13 


B+ K*+ + Lo 


6.0x10"^ 




0.386 


-0.402 


-0.047 


+231 


—11 


B- K*- + /)° 


9.0x10-6 


3.1 


0.109 


-0.335 


0.009 


+0.74 


— O.Od 


B+ K*+ + /)0 


8.5x10-6 




0.109 


-0.335 


0.009 


+2.18 


-0.19 


B- K*- + J/-^ 


3.9x10-3 




0.427 


-0.621 


0.123 






B- D*- + 


2.4x10-2 


-0.03 


0.476 


-0.664 


0.183 








h ^ d transitions: 


As = Ac = 0, A6 = 1 






5- ^ p- + $ 


3.9x10-^ 




0.130 


-0.364 


0.011 






B- K*- + K*^ 


1.4x10-6 


-3.5 


0.112 


-0.352 


0.014 






B~ ^ p- +UJ 


9.1x10-6 


-5.6 


0.083 


-0.298 


0.007 


-0.12 


+0.001 


B+ ^ p+ +UJ 


10.x 10-6 




0.083 


-0.298 


0.007 


+0.18 


-0.002 


B- ^ p- + p^ 


1.3x10-^ 




0.083 


-0.298 


0.007 






B- ^ p- + J/^ 


1.6x10-^ 




0.388 


-0.597 


0.097 






D- — D*^ + ir- 


;i.2xio-'^ 


0.55 


0. 156 


-O.GGO 


0.172 







Table 4 



Observables 


Decay mode 




tree Sz penguin 


pure penguin 


r, 


, CtCp , Cp 


2 22 

f-p ) ^sCfCp , OlgC-j- 


T -f , ai- OLi 


agcj , {asCtCp , asCp) 


(^s^tCp , (o^Cp) , OigC^ 


Pi 


2 2 

C-tCp 1 Cp , C^sC-i 


2 2 2 
Cp , OlgCfCp , OlgC-f. 


A- A 


(agcj , agCtCp , ascl) 


asCtCp , (cisCp) , a^cf 
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